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Abstract 

The reaction — » n + nn is calculated up to order \ 5 ^ 2 i n chiral perturbation theory, where x 
denotes the ratio of the pion to the nucleon mass. Special emphasis is put on the role of nucleon- 
recoil corrections that are the source of contributions with fractional power in \- Using the known 
near threshold production amplitude for jp — > Tt + n as the only input, the total cross section for 
7<i — > n + nn is described very well. A conservative estimate suggests that the theoretical uncertainty 
for the transition operator amounts to 3 % for the computed amplitude near threshold. 



1 Introduction 



Low-energy meson-nucleus reactions are of great interest for they are one of the best tools to deepen 
our understanding of the nuclear few-body problem. In particular, reaction involving the lightest 
member of the Goldstone octet, i.e. the pion, are subject of special experimental and theoretical efforts 
since they can be treated within chiral perturbation theory (ChPT). In this scheme high accuracy 
calculations in connection with reliable error estimates are possible. Indeed, since the pioneering 
works by Weinberg [1] and Gasser and Leutwyler [2] ChPT has developed into a powerful tool for 
investigations of the tttt [3], ttN [4] as well as few nucleon systems [5]. In addition, it was also Weinberg 
who pointed out how to calculate in an equally controlled way pion scattering off as well as inelastic 
reactions on nuclei [6]. There is also a large amount of work in the literature on neutral pion production 
(see Refs. [7, 8] and references therein) as well as on Compton scattering (see Refs. [9] and references 
therein) on the deuteron within ChPT. 

In this paper we will present, for the first time, a complete calculation within ChPT up to order 
^,5/2 £ Qr reac ti on — > TT + nn, where the expansion parameter for the chiral expansion is denoted 
by x = tti^/Mn. Here (Mjy) denotes the pion (nucleon) mass. Data on this channel exist for 
excess energies Q = y/s — (2Mjy + m^) < 20 MeV [10] which allow us to verify our results. The 
calculation is of high theoretical interest, because it provides an important test for our understanding 
of those aspects of ttNN dynamics that are relevant for pion production reactions on the deuteron. 
That understanding is a prerequisite for the reliable extraction of the pion photoproduction amplitude 
on the neutron, commonly done from corresponding deuteron data, but also for the determination of 
the nn scattering length from tt + production data. 

In view of the high accuracy of the data and also because of the high reliability required for the 
extraction of the above mentioned quantities involving neutrons it is now time to critically investigate 
(and avoid whenever possible) the approximations traditionally used in pion reactions on few-nucleon 
systems. Here we will focus on approximations to the pion rescattering contribution as they are 
commonly used in both effective field theory calculations (see Ref. [11] and references therein) as well 
as phenomenological calculations (see Ref. [12] and references therein). Especially for the effective field 
theory calculations one might wonder how recoil corrections could be an issue, for the formalism allows 
for a rigorous expansion in m^/M/v that should build up the recoil corrections perturbatively. However, 
it was stressed recently [13] that the ttNN threshold introduces non-analyticities in the transition 
operators that call for special care: instead of being suppressed by one power in m^/Mjy compared to 
the formally leading rescattering contributions (static term), as one might expect naively, the nucleon 
recoil terms turn out to scale as \Jrn^jMjq relative to the static term as will be demonstrated in this 
paper. The only publication we are aware of, where the recoil corrections were treated properly for 
the near threshold region, is a phenomenological calculation for ^yd — ► ir°d presented in Ref. [14] (so 
far most phenomenological studies concentrated on the A-region, cf. Ref. [15] and references therein). 

Recently the effect of the nucleon recoil on rescattering processes of pions in Trd scattering was 
studied [13, 16, 17] (for previous investigations on the role of the nucleon recoil see Refs. [18]). In 
particular, in Ref. [13] we demonstrated that, at least for the ird system, the nucleon recoil can 
be neglected as long as the two-nucleon intermediate state is Pauli forbidden, while the pion is in 
flight. Thus, in this case the static approximation for the pion exchange is justified. However, as 
soon as the two-nucleon state is Pauli allowed, the nucleon recoil has to be included. In this case it 
turned out that the whole rescattering contribution (i.e. static term + recoil corrections) practically 
canceled completely. It should be stressed that for ird elastic scattering the Pauli allowed two-nucleon 
intermediate states are anyhow suppressed by chiral symmetry. 

In this paper we will study the reaction jd — > ir + nn with special emphasis on the abovementioned 
recoil corrections. In this system the selection rules are such that the S-wave two-nucleon intermediate 
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Figure 1: Diagrams for 7c? — > ir + nn. Shown are one-body terms (diagram (a) and (b) ), as well 
as the corresponding rescattering contribution (c) — all without and with final state interaction. The 
need to include the NN interaction non-perturbatively implies also the inclusion of diagrams with 
the two nucleon pair interacting, while a pion is in flight — this class is shown in diagrams (d). Solid 
straight, wiggly and dashed lines denote nucleons, photons and pions, in order. Filled squares and 
ellipses stand for the various vertices, the hatched area shows the deuteron wave function and the filled 
circle denotes the nn scattering amplitude. Crossed terms (where the external lines are interchanged) 
are not shown explicitly. The tree level yp — > n + n vertex, as it appears in diagram (al) and (a2), 
contributes at leading order (order %°), and order x 1 an d X 2 depending on the one-body operator 
used. Loops start to contribute at order x 2 , the corresponding recoil corrections enter at order x 5 ^ 2 - 



state during the pion rescattering process is allowed by the Pauli principle. In addition, once the 
strength of the one-body operator is fixed to the reaction yp — > tt + u, no free parameters occur in the 
calculation for yd — ► Tr + nn to the order where the leading recoil corrections enter and thus we can 
compare our results to experimental data directly. At the same time we get a better understanding 
of the few-body corrections to yd — ► TT°pn. This reaction will eventually allow one to extract the 
amplitude of yn — ► ir°n complementary to using yd — > ir^d discussed in Ref. [7]. Note that up to now 
no calculation for yd — > Tr + nn exists where the nucleon recoil was properly included. 

Before going into the details some comments are necessary regarding the relevant scales of the 
problem. In the near threshold regime of interest here (excess energies of at most 20 MeV above pion 
production threshold) the outgoing pion momenta are small compared even to the pion mass. Thus, 
in addition to the conventional expansion parameters of ChPT m n /A x and q^/A x , where A x denotes 
the chiral symmetry breaking scale of order of (and often identified with) the nucleon mass, and g 7 
denotes the photon momentum in the center-of-mass system which is of order of the pion mass, we 
can also regard k^/m^ as small, where k n denotes the outgoing pion momentum. In what follows we 
will perform an expansion in two parameters, namely 

Xm = m w /M N and xq = K/m n . 

Obviously, the value of the second parameter depends on the excess energy. The energy regime of 



3 



order 


s-wave 


p-wave 


d-wave 


1 




Am 






i 

X 


1 2 

Xmi Xq 






x 2 


2 2 4 

Xmi XqXitij Xq 

2 5/2 1/2 4 
XmXQ) Xm > Xm Xq 




Xq 



Table 1: Pattern of appearance of the expansion parameters Xm = r^n/^N and xq = k^/m^ on the 
amplitude level for a given order in x f° r the total cross section. The first column shows the order 
parameter whereas the other three columns show the order assignments for the amplitudes of the 
various pion partial waves of relevance. Note, as we here consider the total cross section only, different 
partial waves do not interfere — this was used in the order assignment. 



interest to us are excess energies up to 20 MeV. The maximum value of %q, XQ ax = x/ZQ/mn, possible 
at maximum energy is thus about 1/2. Since this is numerically close to yjxm we use the following 
assignment for the expansion parameter 

X ~ Xm ~ Xq ■ (1) 

In this paper we will concentrate on total cross sections only and thus — up to one important exception — 
XQ appears with even powers only. In table 1 we show the powers of Xm and XQ that appear in the 
amplitude as well as the corresponding order x for the total cross section. The pertinent diagrams 
will be discussed in detail below. Note that the diagrams with irN rescattering (see diagrams (b), 

(c) and (d) in Fig. 1) contribute at order x m as wen as at XmXQ-, Xm" 2 and at xI^Xq- The origin of 
the non-integer power of % are the two-body (irN) and three-body nNN singularities. This issue is 
discussed in detail in section 3. 

The small pion momentum in the exit channel leads to a suppression of higher pion partial waves. 
As can be seen from table 1 we need to consider at most pion d-waves (here and in what follows we 
denote pion partial waves by small letters and NN partial waves by capital letters). On the other 
hand, at excess energies of 20 MeV the maximum two nucleon momentum in units of the pion mass 
is of order 1 and thus there is a priori no suppression of higher NN partial waves. However, since the 
nn phase shifts are only sizable for S- and P- waves at the small energies of relevance, we only include 
the nn final state interaction of those partial waves. 

To summarize the scope of this work, it aims to improve the existing calculations for the reaction 
jd — > ir + nn in the following important aspects: 

• A first complete ChPT calculation for the reaction jd — > n + nn up to order x 5 ^ 2 is presented. 
At this order loops contribute and the non-perturbative character of the NN interaction calls 
for the use of interacting two nucleon Green's functions (leading to up to 3 loop diagrams as 
shown in diagram (d2) of Fig. 1). 

• The leading nucleon recoil is included without approximation. It enters at order x 5 ^ 2 - 
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(a) (b) (c) 

Figure 2: Typical diagrams for nd scattering. Shown are one-body terms (diagram (a) and (b), as 
well as the corresponding rescattering contribution (c)). Crossed terms (where the external pion lines 
are interchanged) are not shown explicitly. 

• As always in effective field theory studies an estimation of the accuracy of the calculation can 
be given. A conservative estimate points at an accuracy of 2 % for the few-body corrections 
to the amplitude near threshold which is of the same order as the uncertainty of the input 
quantity Eq + — the invariant electric dipole amplitude for *yp — ► ir + n at threshold. Adding the 
two uncertainties in quadrature we arrive at a total uncertainty of 3 % for the full transition 
operator. In this work we use phenomenological AW wavefunctions and thus we are not in the 
position to estimate the uncertainty of the complete matrix elements (see the corresponding 
discussion in the summary). 

• For the first time pion p- and d-waves as well as the final state interaction in the NN P-waves 
are included in a calculation for the near threshold regime. Note that the significance of NN 
P-waves even at photon energies below 20 MeV was established long ago [19], however, so far 
they were only considered as plane waves in the impulse term (thus only through diagram (al) 
of Fig. 1). 

Our paper is structured as follows: in the next section we briefly review the central findings of 
Ref. [13]. The three-body dynamics is discussed in detail in section 3. The same formalism as 
described for ird scattering in section 2 is applied to the reaction jd — * ir + nn in section 4. In section 5 
we present the results followed by a brief summary. The details of the calculation for the various 
diagrams are given in appendices. 

2 Remarks on the ird system 

To keep the paper self-contained, in this section we briefly review the findings of Ref. [13]. The 
relevant diagrams for ird scattering are shown in Fig. 2. Diagram (a) denotes the tree level one-body 
contribution, diagram (b) the loop correction to the one-body piece and diagram (c) the rescattering 
term. It should be clear that the Pauli principle calls for a consistent simultaneous treatment of (b) 
and (c), for an interchange of the two nucleons while the pion is in flight (this intermediate state is 
marked by the perpendicular line in Fig. 2) transforms one diagram into the other. 
The calculations of Ref. [13] were based on the effective ttN vertex, 

Kn = 9+ + t g- , (2) 

where a, b, and c are the Cartesian pion indices. We note that this vertex leads, up to corrections of 
0(m n /Mi\i), to identical results for the nd scattering length as the leading order chirally symmetric irN 
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interaction (the Weinberg-Tomozawa term) if we choose g_ = —y/2m n /(2f%) and g + = 0. However, 
to keep the expressions simple we will use the notation of Ref. [13]. 

The main issue of Ref. [13] was to properly isolate the single nucleon contribution (the one that 
would be measured in ttN scattering) from the few body corrections. It is clear that part of diagram (b) 
contributes to the former and part to the latter. As outlined in Ref. [13] the proper prescription to 
separate these two pieces is to add to the tree level scattering (diagram (a)) the single nucleon one- 
loop contribution on the free nucleon at rest — this sum is the expression for the ttN scattering length, 
a (i-body)_ ip^g same i OQ p neec i s to be subtracted from the full contribution depicted in diagram (b). 
This procedure at the same time renders the expression for the loop finite. This difference is a true 
two-nucleon operator. To obtain more symmetric, easier to interpret results we subtracted from the 
rescattering piece the expression for the static exchange, i.e. the contribution from diagram (c) of 
Fig. 2 calculated in the static limit. In order to leave the final result unchanged, this contribution 
needs to be added as d^Q (and fl/vxo) to the expression for the ird scattering length. Thus, the full 
result for the ird scattering length reads 

a = a (i-**) + a Jg + a (rec) + , (3) 

where the individual contributions for the static (st), the recoil (rec) and the NLO corrections to the 
static term are given by 

a { $={9%-9 2 -)h ; a^= g \l ++g ll- ; a^-^(<£-/)/ . (4) 

The integrals denoted by Iq and I±, given explicitly in Ref. [13], are evaluated numerically using the 
deuteron wave functions from the Bonn potential [20]: 

I + = -0.88 I , J_ = -0.19 I . 

These numbers clearly reflect the claim made above: for the ird system the isovector irN interaction, 
proportional to <?_ , leads to a Pauli blocked intermediate state and numerically we find that the recoil 
corrections lead only to a 20% correction (Integral J_ is numerically small, /_ <C Iq.). When the NLO 
pieces mentioned above are added, the correction resulting from recoil and NLO terms together is only 
4% and thus the static pion exchange is a good approximation for the total rescattering contribution. 
On the other hand, the isoscalar irN interaction leads to a Pauli allowed intermediate state and in this 
case the recoil corrections cancel 90% of the static exchange (J+ is large and negative, (Io + I+) <C Io). 
The inclusion of the NLO piece in this case further reduces the total rescattering contribution down 
to 3 % of the static term. Therefore, in this case estimating the total rescattering contribution by the 
static exchange only is a very poor approximation. 



3 The role of the nNN cuts 

In this section we discuss the role of the three-nucleon cuts in more general terms. All arguments pre- 
sented for the three-body dynamics apply to pion-nuclear reactions in general (as, e.g., ird scattering), 
however, in this section we will use only the reaction 7c? — ► ir + nn as illustrative example. 

Typical diagrams that contain a irNN intermediate state (cf., e.g., diagram (c2) of Fig. 1 and the 
corresponding Eq. (33) in appendix A), can be cast into the following form by a proper choice of 
variables (and dropping terms of higher order in iti^/Mn) 

[ dK^KldPRPl f{KlPp 
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2 

Figure 3: Notation used in the discussion of the ttNN intermediate state. The ellipses labeled 0\ 
and O2 denote the transition operators for the various possible reactions; e.g. for 7c? — > ir + nn the 
operators 0\ and O2 refer to 7a! — > 7r + nn and 7r + nn — > 7r + nn, respectively. 



where iT^ denotes the pion momentum and Pr denotes the relative momentum of the two nucleons 
while the pion is in flight (c.f. Fig. 3). Using the static approximation means to expand the denom- 
inator of the integral in Eq. (5) in powers of m^jM^ before evaluation of the integral. Thus, we get 
in the static approximation* 1 

Astatic) (n , _ f dK«KldP R Pl ml Pi) (rn^A 

KnN {Q) ~J (2vr)6 Q-K$/(2m„)+iO + U {M N ) ' W 

We will demonstrate now analytically (and in section 5 numerically) that this procedure misses an 
important contribution to I n NN- 

Above the pion production threshold the denominator in Eq. (5) has a three-body singularity that 
leads to an imaginary part. This imaginary part of I^nn can be calculated by the replacement 

(Q - Kl/{2m«) - P 2 R /M N + iO)- 1 - -mS{Q - ^/(2m 7r ) - P 2 R /M N ) . 

Thus we find 

&(Q) = ~™™*J ^^K^AQ - P 2 R /M N ),P 2 R )^2m n {Q - P 2 R /M N ) . (7) 

From this formula we see that the imaginary parts of the pion loops in diagram (b2) and (c2) of 
Fig. 1 are finite and lead to a strongly energy-dependent contribution. As long as the momentum- 
dependence of the function / can be neglected, this part of the amplitude grows like Q 2 , i.e. like the 
three-body phase space. These are the contributions at order X^Xq given in table 1. Note that 
the corresponding amplitudes without final state interaction ((bl) and (cl) of Fig. 1) do not have a 
three-body cut but only a two-body singularity. Thus their imaginary parts scale as XmXQ- 
On the other hand, the unitarity cut contribution of 1^^°^ reads 

& StattC \Q) = -i^y/^Q J ^ff(2rn«Q,P 2 R ) ■ 

The most remarkable difference between and I^^ taU< ^ is that in contrast to the former the 

latter scales as y/Q, i.e. like the two body phase space for all diagrams (even for those with final state 
interaction) and therefore shows an energy dependence that is completely wrong — and at variance 
even with perturbative three-body unitarity. 

#1 Note, we here use the phrase 'static approximation' in a quite broad sense in that we also allow for the inclusion of 
correction terms analytic in \- In phenomenological studies those corrections are normally dropped. 
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For values of Pr with Q > P R /M^ the integral I^^(Q) (cf. Eq. (7)) contributes to the imaginary 
part of IjrNNiQ) (cf. Eq. (5)). To evaluate the contribution of the cut to the real part, the square 
root needs to be analytically continued to negative values of its argument through 

y/2m n (Q - P 2 R /M N ) - i^2m^(P 2 R /M N - Q) . 
To demonstrate explicitly the impact of this, let us consider the case Q = 0. Then we get 

iSi(O) = «m*^ J d -^f{-2m,Pl/M N ,Pl)P R . (8) 

The corresponding expression for I^ t ^ tatlc \d) vanishes. On the other hand, the expression for the 
leading static approximation (cf. Eq. (6)) gives at threshold 

(0) = "2m. f ^^ ml Pi) • (9) 

The static approximation Eq. (9) only acquires corrections analytic in (m^/M/v) and thus misses the 
contribution of Eq. (8). In fact, Eq. (8) corresponds to the threshold value of the mentioned non- 
analytic contribution from the three-body intermediate state that is dropped in the static approxi- 
mation — the extension to arbitrary values of Q is straightforward. However, the contribution from 
^TrNN ^ s significant as one can see from a naive dimensional analysis where all momenta — even those 
in the integral measure — are replaced by their typical values: 

T (cut) j m n (static) / 1fl s 

nNN y A^/y" ' ^ ' 

as claimed in the introduction. Therefore, in general the static approximation is to be avoided! As 
one can see, the contributions from the nucleon recoil through the three-body singularities to both 
the real and the imaginary part of the amplitude appear to be down by y/x compared to the leading 

loop contribution (the static piece). To be more explicit they contribute at orders Xm 2 and Xm 2 XQ 
(cf. table 1). On the other hand, if we expand the propagator in Eq. (5) before the integration, we 
only get terms analytic in the pion mass (order x corrections to the static piece) and miss the most 
prominent correction. 

Now we are in the position to discuss in more detail the conjecture presented in Ref. [13] as well as 
in the previous section, namely that in general for all those diagrams, where the S-wave two-nucleon 
intermediate state that appears while the pion is in flight is forbidden by the Pauli principle, the 
various recoil corrections largely cancel, while in case of Pauli allowed S-wave intermediate states 
they add coherently. In the latter case the net effect of the rescattering diagrams was claimed to 
be small due to a destructive interference between the recoil corrections and the leading rescattering 
contribution (the static term — cf. Eq. (6)). We will now discuss the two cases in the light of the 
discussion above. The essential observation is that the recoil corrections are the analytic continuation 
of the imaginary parts related to on-shell itNN intermediate states. 

Pauli forbidden intermediate states: in this case the imaginary contributions stemming from the 
three-body unitarity cut in diagrams of the type of (b2) and those from the type of (c2) of Fig. 1 
need to cancel exactly, for the corresponding irNN state is forbidden. As a consequence, there will 
also be a cancellation for the analytic continuations (see Eq. (7)), i.e. for the corresponding principal 
value (PV) integrals, and thus the recoil corrections necessarily cancel to a large degree. Since this 
statement is based solely on the Pauli principle it must hold for all reactions. 
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Figure 4: Diagrams that contribute to the one-body production operator in diagrams (al) and (a2) 
of Fig. 1. Only tree level diagrams that contribute to 7p — > 7r + n up to order \ 2 are shown. The dots 
denote loops that contibute to the s-wave amplitude at oder \ 2 as those shown in Fig. 5 (b) and (c). 



Pauli allowed intermediate states: contrary to the first option in this case the corresponding uni- 
tarity cut parts of the diagrams (b2) and (c2) of Fig. 1 add coherently and, as above, the same is to 
be true for their analytic continuation. It is not possible to claim in general that the recoil contribu- 
tion from the PV integrals completely cancels the whole static term. However, as illustrated by the 
estimate of Eq. (10), the recoil corrections tend to be of the order of magnitude of the static term and 
tend to interfere destructively with it. Here we refer to the discussion in section 5 and to Fig. 10. 

As we discussed in the previous section, the leading operator for ttN scattering acting on a deuteron 
leads to a irNN state, where the S'-wave is forbidden for the NN-p&ii by the Pauli principle. On 
the other hand, as will be outlined below, the leading operator for — > irN, when acting on a 
deuteron field, leads to a irNN state, where the AW-pair is allowed to be in an S-wave. Thus, in 
the reactions ird — ► ird, ird — ► jNN, and jd — > 7r°d the static approximation indeed accounts for a 
significant fraction of the few-body corrections, whereas in 7a! — > ir + nn and 7a! — > 7c?, when evaluated 
near the pion threshold, the static approximation is expected to work quite poorly. References for the 
corresponding calculations are given in the introduction. 

4 The reaction "yd — » 7r + nn 

We will start with a discussion of the one-body contributions to be used in the diagrams of type (al) 
and (a2) of Fig. 1. 

At threshold the ir'yNN vertex at leading order (Xm) an< ^ a * next-to-leading order (Xm) is given 
by the so-called Kroll-Ruderman (KR) term [21] and its recoil correction # 2 , 

V£U = ieg nN (l - (6 7 • a) e 3ab r b , (11) 

where e* 7 denotes the photon polarization and uj w is the energy of the outgoing pion. The corresponding 
diagram is shown as diagram (a) in Fig. 4. Note that we use g n N = 13.4 in the calculation and the 
charge e is normalized such that the fine structure constant is given as e 2 /(4-7r) = a = 1/137. This 
vertex contributes to the one-nucleon operator (impulse term) as shown in diagram (al) and (a2) of 
Fig. 1 and also provides the production vertex for the virtual pion to be rescattered. As the ttjNN 
vertex is both spin and isospin dependent, now, in contrast to ird scattering, for the rescattering 
contributions the two-nucleon intermediate state that occurs while the pion is in flight is allowed by 
the Pauli principle to be in an S'-wave. 

* 2 A list of the vertices relevant for chiral perturbation theory calculations can be found in Ref. [22]; however, in contrast 
to that reference we use the spinor normalization uu — 2Mn ■ In addition we used the Goldberger-Treiman relation to 
replace pa//tt by p^jv/Mjv. 
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(a) 



(b) 



(c) 



Figure 5: Typical pion loops that contribute to 7p — > 7r + n. Crossed terms for diagrams (a) and (b) 
are not shown explicitly. 



According to table 1 at order \ we a l so need to consider the leading correction in Xq to the s— 
wave as well as the leading p-wave contribution which is of order of XQ- Both are provided by the 
diagram where the photon couples to the pion in flight corresponding to diagram (b) of Fig. 4. The 
corresponding expression for the effective ttjNN vertex reads 

V^NN = ~^nB ■ (K - ? 7 ) J]' K r T e 3ab r b , (12) 

where uj n denotes the energy of the outgoing pion. This vertex contributes to the leading p-waves 
(order xq for the amplitude — order Xq = X f° r the total cross section), the leading energy dependence 
of the s-wave (order Xq f° r the amplitude as well as the cross section, for this term can interfere with 
the contribution of order x°), as well as the leading ci-waves (order Xq f° r the amplitude and thus 
order Xq = X 2 f° r the cross section). 

At order Xm Pi° n loops start to contribute to the s-wave part of the ^yp — ► 7r + n amplitude (see 
Fig. 5) and thus to the reaction jd — > ir + nn as well. The only additional vertex needed for the 
evaluation of these loops is given in Eq. (2) with = - v / 2m 7r /(2/2) and g + = 0. The imaginary 
part of diagram (a) in Fig. 5 contributes at XmXQ whereas the imaginary part of diagram (b) starts 
to contribute only at higher order (XmXg) an d is thus not explicitly included in the calculation. In 
contrast to the imaginary parts the real parts of some of those pion loops on a single nucleon are 
divergent and need to be regularized (e.g. diagram (c) of Fig. 5) # 3 . In the calculation of the reaction 
'yd — > -K + nn we use a prescription similar to that used for ird scattering, as described in the previous 
section. In practice this means to replace the expression of Eq. (11) by 

V eff = iKE 0+ (e~ r a)e 3ab T b , (13) 

where k = 4tt\/~2 (m^ + M^). The experimental value of Eq + from the reaction 7^ — ► ir + n is [25]: 

E 0+ = (28.06 ± 0.27 ± 0.45) x lO^m" 1 . 

This value coincides within errors with the result of ChPT up to order xfn [26]: 

E 0+ = (28.2 ± 0.6) x lO^m" 1 . (14) 

Although this calculation was performed to higher order than what we aim at here, we will use the 
latter value as input quantity for the threshold value of Eq + since the contribution of order Xm are °f 
the order of the assigned uncertainty. 

#3 For a complete discussion of all relevant loops we refer to Refs. [23, 24]. 
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It is straightforward to show that the s-wave contribution derived from Eq. (12) has the same 
operator structure as Eq. (13) and we may simply include its effect by replacing Eq + in Eq. (13) by 
E 0+ (k 2 ) defined up to order Xq through 

E 0+ (kl) = E 0+ + E> + ^ + E^ . (15) 

The leading contribution to both coefficients Eq + and Eq + can be calculated directly from Eq. (12). 
At next to leading order the energy dependence of the ^ttNN vertex enters as well (c.f. Eq. (11)). 
Especially, we find 



^ = -^(i + ^)=-mxiir«™.-' 



where |(f 7 | = m n (l — m n /(8MN) + 0{\ 2 )) was evaluated in the center of mass system of the reaction 
'yd — > iT + nn at threshold. The slope parameter Eq + is calculated up to its leading and next-to-leading 
order. However, we anyhow assign an uncertainty of 10 % to this quantity, for higher order corrections 
are known to be enhanced by the A resonance [27]. The given value for E^ + compares well to that 
from the dispersive analysis of the Mainz group [28] and is consistent with that used as input in 
Ref. [25] to extract E$ + from the data. 

At order x 2 f° r ° ur reaction, additional diagrams contribute to the transition 'yp — > ir + n where the 
photon gets absorbed on the nucleon, e.g. through the magnetic moment, followed by pion emission 
(diagrams (c) in Fig. 4). Thus, they are accompanied by a one-nucleon intermediate state. These 
diagrams, when used as the one-body operator in diagram (a2) of Fig. 1, acquire a two-nucleon cut 
that leads to an imaginary amplitude even at threshold due to the kinematically allowed transition 
■yd — > pn followed by pn — > ir + nn. However, this two-nucleon cut introduces a new (large) momentum 
scale p ~ ^m^Mjy into the problem that calls for special care. According to the counting rules for pion 
production in NN collisions [29], the contribution of this two-nucleon cut is therefore suppressed by 
Xm 2 compared to the rescattering diagram (c2) of Fig. 1. It is therefore justified within the ordering 
scheme used to replace the two-nucleon propagator by its static limit. The resulting contribution 
from the magnetic couplings to the s-wave pion production on the nucleon is already included in the 
effective operator of Eq. (13). For the pion p-waves we get the following expression for the sum of the 
s- and u-channel contributions of Fig. 4(c) 



y(c) _ ieg n N 



^ NN " V2M n m T 



+ , (16) 



where 



A = i{n p - ii n )K ■ (e* 7 x <f 7 ) , (17) 
B = -2 (e 7 • p) K + (fip + fin) ((e 7 • k n )q 7 - (<f T • K)e^j . (18) 

Here p denotes the momentum of the incoming nucleon and fi p = 2.79 and fi n = —1.91 denote the 
magnetic moments of the proton and the neutron, respectively. The operator V^yNN & Yen by Eq. (16) 
was evaluated directly for the channel 'yp — > 7r + n. As a consequence the expression does not show any 
explicit isospin dependence and an isospin factor of y/2 appears. 

In Ref. [27] it is shown that the pion p-waves converge quite slowly: the next-to-leading order 
correction to the dominant p-wave multipoles M]+ and M x - change the leading result by a factor of 
2. The reason for this sizable correction is the large numerical value of {fi p — fi n ) = 4.7. On the other 
hand, the contribution of the A-tevm to 'yd — > -K + nn is suppressed. One reason can be read off from 
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Figure 6: Diagrams where the photon couples to the structure of the deuteron. These diagrams start 
to contribute only at order x 3 and are therefore not considered here. 



Eq. (16) almost directly: since the ^4-term is a scalar in spin space, it will not change the total spin 
of the two nucleon system, when used as the one-body operator in diagram (a) of Fig. 1. Thus the 
final nn state is in a spin-triplet state. However, in the near threshold regime of interest here, the 
simultaneous appearance of an nn P-wave (nn spin-triplet states are to have odd angular momenta 
as a consequence of the Pauli principle) and a pion p-wave is suppressed. In addition, it turns out that 
for the total cross section the A-term does not interfere with the leading pion p-wave contributions 
(the term proportional to <f 7 in Eq. (12)) which leads to an additional suppression. For details on the 
latter point we refer to the explicit expressions given in appendix A. One should also note that like 
for the slope of the s-wave amplitude, the A-isobar gives potentially large corrections to the p-wave 
amplitudes [27]. Thus we assign an uncertainty of 10% also to those. 

At order \ 2 , pion rescattering diagrams (see diagrams (b), (c) and (d) in Fig. 1) start to contribute. 
Here the same diagrams contribute as previously discussed for ird scattering, but with the first irN 
interaction being replaced by the — > irN transition vertex. These diagrams are depicted in Fig. 1 
(b2) and (c2). However, there is in addition a whole new class of diagrams, namely those without final 
state interaction (depicted as (bl) and (cl)). Another important difference to nd scattering is that, as 
already mentioned above, the S-wave two-nucleon state — while the pion is in flight — is now allowed 
by the Pauli principle. This has two consequences: first of all we expect that the static approximation 
will not work for the rescattering contribution and, as a second consequence, now nothing prevents the 
two-nucleon system to interact while the pion is in flight. As the NN interaction is to be taken into 
account to all orders, these diagrams are also potentially important. The latter statement becomes 
especially clear when observing that the scattering length in the nn channel is quite large; note that 
in the limit of an infinite scattering length the diagrams (dl) and (d2) acquire a logarithmic infrared 
divergence. In fact, then the two-nucleon propagator behaves as if it would describe the propagation 
of a massless particle* 4 . 

Starting from the vertices given by Eq. (11) and in the appendix of Ref. [22] it is straightforward 
to write down the corresponding matrix elements for the diagrams shown in Fig. 1. Note that all 
diagrams are evaluated in Coulomb gauge. This is a standard choice in ChPT calculations, for many 
diagrams with single photons are relegated to higher orders, such as those where the photon couples 
to the charge of the nucleon and of the deuteron, respectively. 

Diagrams, where the photon couples to a rescattered pion in flight turn out to be strongly suppressed 
numerically compared to those shown in Fig. 1. The reason for this is twofold: their imaginary part 
is suppressed as a consequence of gauge invariance that forces the photon-pion coupling to be of the 
order of the (small) pion momentum in the loop and the appearance of a second pion propagator 
leads to an additional suppression. This is in complete analogy to nd scattering as discussed in detail 
in Ref. [31]. In addition gauge invariance in principle also calls for the inclusion of diagrams, where 

* 4 In a different scheme and for a different reaction this behavior was studied in Ref. [30] . 
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Figure 7: Total cross section of the reaction -yd — ► ■n + nn at LO (dashed line), NLO (dash-dotted line) 
and x 5 / 2 order (solid line) together with experimental data from Ref. [10]. 



the photon couples to the internal structure of the deuteron. Typical representatives of this class are 
shown in figure 6. However, it is easy to show that these diagrams are suppressed by at least one 
power of x compared to the leading rescattering contribution and therefore start to contribute at most 
at order x 3 to 'yd — > 7r + nn* 5 . Therefore these diagrams are not included in the calculation. As a 
consequence, the diagrams given in Fig. 1 (together with Figs. 4 and 5) are all that contribute to the 
given order. 

Details on how the various diagrams are evaluated are given in the appendix A. Especially, there it 
is explained how we dealt with the three-body singularities that occur in diagrams (c2), (b2), (dl), and 
(d2). We want to mention that we evaluated the loop diagrams for non-relativistic pions (see appendix 
A) for that largely simplified the numerics. We checked by direct evaluation of diagrams (cl) and (c2) 
for threshold kinematics that switching to relativistic pions changes the individual contributions to 
the amplitude by less than 4 %. 

5 Results and discussion 

In our calculations we use the standard values for the various constants, namely /„- = 92.4 MeV, 
m n = 139.57 MeV (only the charged pions contribute to the order we are working) and Mn = 938.27 
MeV. The deuteron wave function and also the scattering wave functions needed for the nn scattering 
amplitudes (in the 1 S'o, 3 Po, 3 -Pi, and 3 P2 partial waves) are generated from the (charge-dependent) 

* 5 For this estimate we used m n as typical momentum in the deuteron. A more accurate estimate of this quantity would 
be 7 = \/EbMn- If we were to use this for the power counting it would yield a suppression of the order of 7/Mjv ~ 0.3x- 
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AE y [MeV] 

Figure 8: Importance of the slope parameter Eq + denned in Eq. (15): The solid line shows the 
full result, whereas the long dashed line is produced with Eq + = 0. The two dotted lines show the 
estimated uncertainty of full result as described in the text. The data are from Ref. [10]. 



CD-Bonn potential [32]. Specifically, for the former we employ the analytical representation of the 
deuteron wave function provided in that reference because it allows us to perform some integrations 
analytically — see appendix B. For the same reason the scattering wave functions are computed from 
rank-1 separable representations of the CD-Bonn model, constructed along the lines of Ref. [33] uti- 
lizing the so-called Ernst-Shakin-Thaler method [34], cf. also appendix B. Note that the scattering 
length predicted by the CD-Bonn potential for the 1 So partial wave in the nn system is a nn = —18.97 
fm [32], which is in line with most of the recent experimental information [35, 36] (note, the analy- 
sis of Ref. [37] gave a significantly lower value). We want to point out that in our calculation the 
contribution of the deuteron D-wave is included. This contribution to the leading diagrams is rather 
important, because it guarantees the correct normalization of the S-wave component for the potential 
used. 

Although it would be desirable to use wave functions consistent with the transition operator, as 
provided, e.g., in Ref. [38], various studies show a low sensitivity to the wave function employed (see, 
e.g., [31]). In addition, the current operators have not yet been worked out to the same order and 
within the same unitary transformation as the NN wave functions. We postpone this to a future 
study. 

The reaction jd — > iT + nn was calculated in the DWBA approximation by many authors in the 
middle of the seventies (cf., e.g., the review paper by Lagct [19] and references therein). The approach 
used in those investigations corresponds to the evaluation of diagrams (al) and (a2) of Fig. 1 using 
the jp — * ir + n vertex of Eq. (11); thus, in our language those were incomplete calculations up to 
next-to-leading order (x), because the energy dependence of the r ynNN vertex and higher pion 
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AE y [MeV] 

Figure 9: The role of the higher pion partial waves and effect of the NN interaction in the P-wave. 
The dashed line shows the calculation where only the S-wave interactions are taken into account, 
the dotted line includes also the higher pion partial waves and the solid line includes in addition the 
contribution of the NN final state interaction in the three possible included in diagram (a2) of Fig. 1. 



partial waves were neglected. Since the most important contribution to the jp — > 7r + n operator in 
the near threshold region originates from the Kroll-Rudermann operator and its first correction (see 
Eq. (11)), which is known, and the convergence of diagram (a2) of Fig. 1 is provided by the universal 
fall-off of the deuteron wave function for small momenta (fixed by the deuteron binding energy) all 
those calculations led to very similar results (that we reproduce). To improve the calculations, some 
authors used instead of the prefactors of Eq. (11) the experimental input for Eo + (k^) at threshold 
as the strength parameter for the one-body term (see, e.g., Ref. [10]). This, however, corresponds to 
an incomplete order-x 2 calculation, as discussed above. In addition, none of the works reported in 
Ref. [19] considered the D-wave of the deuteron consistently. The slope of -Eo+(& 2 )) pion p-waves, 
or the NN final state interaction in the P-waves were also not considered. The only attempt to 
improve the mentioned calculations via inclusion of a pion rescattering contribution was made in 
Ref. [39], where diagram (c2) of Fig. 1 was evaluated in the static limit, i.e. without nucleon recoil. 
This contribution was found to be large, amounting to an increase of around 10% of the total cross 
section (see the discussion of this question in Ref. [19]). However, as we stressed above, the static 
approximation is expected to work very poorly in this reaction. 

Our results at order x°; X an d combined order x 2 an d X 5 ^ 2 are shown in Fig. 7 by the dashed, 
dash-dotted and solid lines, respectively. The data are taken from Ref. [10] and the energy is measured 
in terms of AP 7 = (l+m n /(2Mj\[))Q — the photon lab energy subtracted by its threshold value. Up 
to order x 1 the result is parameter free, i.e. it is determined by the value of g^N only. For the value of 
Po+(fc 2 -) at threshold that is needed as the only input quantity at order x 2 , we use the central values 
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Figure 10: Quality of the static approximation for diagram (c2) of Fig. 1: The dash-dotted line shows 
the full calculation for diagram (c2), the dashed one that for diagram (b2), and the solid line the sum 
of the two. On the other hand, the result for diagram (c2) evaluated in the static approximation is 
given by the dotted line. 



of the ChPT calculation of Ref. [26] (that agrees with the experimental number, as mentioned in the 
previous section). 

To illustrate the importance of the slope Pq + we show in Fig. 8 both the full result as well as the 
one we get for a vanishing slope. For AP 7 larger than 10 MeV we observe a quite large sensitivity to 
the energy dependence of Eo + (k%). 

The counting scheme of Eq. (1) calls for an inclusion of p— and d-waves for the pion. The explicit 
expressions for the vertices are given in Eqs. (12) and (16). These vertices are to be considered 
embedded in diagrams (al) and (a2) for Fig. 1 only. The contribution of the higher pion partial waves 
can be read off the difference of the dashed and the dotted line in Fig. 9. 

The significance of NN P-waves, even for AP 7 < 20 MeV, was established long ago [19], however, 
so far they were only included as plane waves (through diagram (al) of Fig. 1). In this work for 
the first time the final state interaction in the NN P-waves is considered as well: we include the 
NN interaction in the three possible P-waves ( 3 Po, 3 Pi and 3 P2) in diagram (a2). The effect of the 
nucleon P-waves on the total cross section is shown in Fig. 9. It is interesting to note that although 
the resulting effect of the P-wave NN interaction is not significant for the total cross section the 
individual contributions of different partial waves are quite large (e.g., the 3 Pi partial wave enhances 
the cross section by about 15% at AP 7 = 20 MeV). It turns out that the contribution of the 3 Pi 
partial wave, which is enhanced due to constructive interference of the diagram (al) and (a2) of Fig. 1 
for this particular NN partial wave, cancels to a large extent the contributions of the 3 Po and 3 P2 
partial waves. Note that the interference pattern of the various NN P-waves can be traced to the 
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fact that 3 Pi is repulsive, whereas 3 P an d 3 P2 are attractive for small relative momenta. We have 
to mention that in the present calculation we neglect the coupling of the 3 P2 partial wave to the 3 i 7 2- 
However, in view of the cancellation mentioned above this coupling might play a role, specifically 
because the mixing parameter €2 is already about one third of the corresponding 3 P2 phase shift at 
NN energies corresponding to Q = 20 MeV. Also the NN 1 Z?2 partial wave, which is likewise not 
included in the present calculation, should start to play a role with increasing energy. As mentioned 
above, the NN wavefunctions used are not consistent with the transition operator. As long as this is 
not the case we do not think that including these probably minor effects is worth the effort. However, 
once wave functions consistent with the chiral counting are used, also the mentioned approximations 
need to be abandoned. 

To summarize the results, we have calculated higher order chiral corrections to the impulse term 
used, e.g., in Ref. [19]. We found that due to a compensation amongst the effects of the various NN 
P-waves, their net effect on the total cross section is negligible. In addition, the effect of the energy 
dependence of the pion s-wave production (parametrized by the slope E' 0+ ) and of the pion p-waves 
compensated each other largerly as well. In view of this findings, the good agreement of the results 
of, e.g., Ref. [19] for the total cross section for photon energies above 10 MeV with the data should be 
considered accidental. However, for differential observables and especially for polarization observables 
we expect sizable effects from the chiral corrections calculated in this work. 

Let us now discuss some of the results in detail. We found that all rescattering contributions of 
order x 2 an d of order x 5 / 2 contribute with similar strength. This is true also for those diagrams 
where the NN interaction appears inside the loop (diagrams (dl) and (d2) of Fig. 1). This clearly 
demonstrates the need to take into account the NN S-wave interaction non-perturbatively also in 
intermediate states (as it is done routinely in three-nucleon calculations anyway). 

Let us now compare in detail the results of different diagrams at the threshold. Here contributions 
from higher partial waves as well as those from the slope parameter E' 0+ vanish. Therefore the uncer- 
tainty of the calculation in this regime can be estimated much more accurately than at higher energies. 
The results of the calculation of all possible contributions to the reaction amplitude at threshold are 
given in table 2. Evidently, the total contribution from few-body (rescattering) corrections to the 
full amplitude at threshold is about 5 % of the contribution from the distorted wave impulse ap- 
proximation. One can also see from this table that the net contribution of the orders Xm an d Xm, 2 
cancels largely the NLO contribution (see also the corresponding results in Fig. 7 at low energies). 
Due to this one might speculate that corrections to this result from higher orders, especially from 
X%ii may influence the results stronger than suggested by the power counting. However we think that 
the uncertainty of the calculation at low energies is indeed within the uncertainty for Eq + given by 
Eq. (14). Actually, no additional diagrams to those shown in Fig. 1 appear at order Xm- Thus, the 
contributions to the amplitude at this order originate basically from three sources: 

1. NLO correction to the Kroll-Rudermann vertex when being used as a one-body operator in 
diagrams (b),(c) and (d) of Fig. 1 and the recoil corrections to the Weinberg-Tomozawa term. 
These corrections are of order (m n /2MN) • 5% = 0.4%. 

2. Relativistic corrections to the diagrams with pion rescattering. For diagram (c2) the relativistic 
correction gives only about 4% (to the 3 % contribution to the amplitude from diagram (c2)). 
Since there are no reasons a priori which could enhance the corresponding corrections to the 
other rescattering diagrams, we assign an uncertainty of order of 0.5 % to this effect. 

3. The contributions from the coupling to the deuteron structure and the nn scattering T-matrix. 
These are the hardest to estimate. However, as argued above, they are expected to be numerically 
quite small (see footnote # 6). 
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operator 


order 


diagrams 


contribution 




Y° 
Am 


(al) + (a2) 


1 


one-body 


V 1 

Am 


(al) + (a2) 

V / V / 


-0.07 




V 2 
Am 


(al)+(a2j 


+0.028 




5/2 
Am 


(b2) 


-0.016 




a/ 2 ^/ 5/2 
Am'A" 1 


(cl) + (c2) 


+U.Uoy 


few-body 


XmiXrn 








(dl) 


+0.008 




2 5/2 


(d2) 


+0.024 



Table 2: Contributions of different diagrams and operators to the reaction amplitude at threshold are 
given at different orders. The results in table are normalized to that of the calculation of diagrams (al) 
+ (a2) with the leading order Kroll-Rudermann vertex. Note that diagram (bl) does not contribute 
at threshold. 



To summarize, we assign an uncertainty of 2 % to the few-body corrections of the amplitude. On the 
other hand the contribution at order Xm t° ^0+ at threshold is also about 2 % [26]. Adding these 
two uncertainties in quadrature we end up with a total uncertainty of 3 % for the transition operator 
near threshold. Unfortunately we are not in the position to estimate the uncertainty from the NN 
interaction used — this has to be postponed until a calculation with fully consistent wave functions is 
performed. 

In Fig. 8 the resulting uncertainty for the full calculation is shown by the two dotted lines. This 
range of uncertainties contains besides the 3 % just discussed also the 10 % uncertainties on both the 
pion p-waves as well as the slope parameter. 

To show that the static exchange is indeed a poor approximation to the exact result for diagram (c2) 
of Fig. 1, as conjectured above, we compare in Fig. 10 the results for the static approximation with 
those from the exact calculation as well as with those for the one-body term (b2). As can be clearly 
seen, the static approximation fails to describe the full result in both strength as well as energy 
dependence. Evidently, for the exact calculation the total contribution of the sum of the results 
for (b2) and (c2) is extremely small near threshold, i.e. in the region where the real parts of both 
diagrams dominate. With increasing energy the contribution from the sum of these two diagrams 
increases rapidly. The reason for this effect is that at higher energies the role of the imaginary parts 
of both diagrams, which contribute coherently for the Pauli allowed AW states, is growing rapidly. 
The resulting contribution for the sum of both diagrams of Figs. 2 (b2) and 2 (c2) to the total cross 
section does not exceed 4%. 
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6 Summary 



In this paper we presented for the first time a ChPT calculation for the reaction 7c? — ► ir + nn. We 
calculated the diagrams displayed in Fig. 1 — keeping explicitly the nucleon-recoil for the intermediate 
states. This corresponds to a complete calculation up to order x 5 ^ 2 > where x = vn^jM^. 

The results of the full calculation are shown in Fig. 7 by the solid line. A good agreement between 
theory and experiment at low energies is obtained without any free parameter. The only input 
parameter was the threshold value for Eq+, taken from an N 3 LO calculation for 7p — > ir + n [26]. 
Estimated conservatively, the uncertainty for the transition operator of our calculation is about 3 % 
in the amplitude for photon energies AE^ below 5 MeV. Unfortunately we are not in the position to 
estimate the uncertainty from the NN interaction used — this has to be postponed until a calculation 
with fully consistent wave functions is performed. 

We found a strong suppression of the pionic rescattering contributions in comparison to the calcula- 
tion in the frozen-nucleon approximation, i.e. in the static limit. This confirms the suggestion made in 
Ref. [19] and is in line with the general remark of our recent paper [13] where it was conjectured that 
the static limit is not adequate for pion rescattering processes with Pauli-allowed S-wave intermediate 
NN states. 

As a next step a fully consistent chiral perturbation theory calculation for the given reaction should 
be performed. It should include not only the use of wave functions constructed within this framework, 
as provided in Ref. [38], but also a re-evaluation of the transition operator within the same unitary 
transformation used to calculate the wave functions. Only then one can reliably asign a theoretical 
uncertainty to the full calculation and address questions raised in Ref. [42] for the consistency of the 
counting scheme and the scaling of four-nucleon operators. Although we do not expect the latter 
effects to be numerically significant for the reaction considered here, they are potentially important 
for 7c? — > TT°pn, where the leading contribution vanishes. 

Once a fully consistent ChPT calculation and better data are available for the reaction jd — ► -n + nn 
one can consider to extract the energy dependence of Eq + from this reaction, for the total cross section 
is very sensitive to the slope of Eq + , as illustrated in Fig. 8. However, before more definite conclusions 
can be drawn on this issue better data on the total cross section as well as the differential cross section 
in the energy regime of 15-25 MeV is needed. 

The work presented allows one, amongst other things, to address the theoretical uncertainty of 
the nn scattering length extracted from jd — ► ir + nn analogously to the studies of Ref. [11] for 
ir~d — > jnn. These studies are of high interest in the light of the significant differences in the values 
for a nn extracted from different groups in different reactions — c.f., e.g., Refs. [35, 36] and Ref. [37]. 
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Figure 11: Definition of the various vectors of the problem. 



A Matrix elements and observables 

The differential cross section is related to the reaction amplitude via [40] 

da = (2n)W\ Pf - P £l\ Mfi \* ,/K ,/» Jl k * , (19) 
y ' v 1 ; 24g 7 ^/i' /l1 (27r) 3 2Ei (2irf2E 2 (2tt) 3 2u 7T ' v ; 

where the indices 1 and 2 label momenta and energies of the final nucleons, and k w (ui n ) denotes the 
pion three-momentum (energy), and Mfi stands for the sum of all amplitudes given below: Mfi = 
M a i + M a 2 + Mbi + . . . . The photon momentum is denoted by g 7 and \/i = 2Mjy + m n + Q where Mjv 
(%) and Q are the nucleon (pion) mass and the excess energy, respectively. The choice of variables 
is illustrated in Fig. 11. The factor of 1/2 accounts for the two identical nucleons in the final state. 

Let us now present the explicit expressions for the amplitudes corresponding to the diagrams in 
Fig. 1. All calculations are done in the center-of-mass system of the reaction. Note that in the 
calculation of the leading diagrams (al) and (a2) in Fig. 1 the D-wave component of the deuteron 
wave function was taken into account, whereas in all other diagrams which are already suppressed by 
Xm as compared to the leading ones, the inclusion of the D-wave is not necessary. This is also true for 
the diagrams with NN P-wave final state interaction. The loops are evaluated using nonrelativistic 
kinematics even for the pion. 

The calculation of the diagrams of Fig. 1 were done using the standard Feynman rules with deuteron 
vertices (for details see, e.g., the appendix in Ref.[41] ). The explicit expressions for each individual 
term are: 



1. Diagram al. 

M a i = C NNLO xl(<?e 7 ) 

where Cnnlo = 167T\/M/v (Mn + rn n ) Eq + , q± t 2 = pi t 2+Q'y/2, and e 7 and are polarization vectors 
of the photon and the deuteron, respectively. The notation —(1 <-> 2) points that one has to subtract 
the same term with spin and momentum variables of the two nucleons interchanged, in order to get 
properly anti-symmetrized amplitude. 



u(q 2 )(ae D )- 



w(q2) 
' V2 



6 s K ae D, 

12 



<?2 

V2 



X2 



(1 <- 2) , (20) 
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Here u(p) and w(p) are the S-wave and D-wave deuteron wave functions defined in appendix B. 
The expression x\ O x\ corresponds to the spin structure of the final NN pair: 

Xi (^)(^b)-^= X*2 = (eyeb) (xi ^| X2) +^(e 7 x e D ) (x\ x£j (21) 

and analogously for the D-wave part. Here the first and second terms on the r.h.s. are the spin-singlet 
and spin-triplet contributions, respectively (x^X = !)■ 



2. Diagram (a2). 

M a2 = -C NNLO ^rj-r- / 7^-73 < — r u(p+g 7 /2) e 7 e s ) Xi -7= X 2 

. M 1 NN (p+k 7V /2,p 12 ,E 12 ) a p _ ( t o- 2 A\ 

+ a- (r+w^+a " (pW2) !(£i x £d)s l x ' "°T2 X2 J / <22) 

where p\ 2 = (1/2) (pi — p 2 ) and the amplitudes -M^ N are related to the NN partial wave amplitudes 
with angular momentum L as given below in this subsection. The energy of the NN system is 
E12 = p\ 2 /Mn ■ The formula for M a 2 as given shows the contribution from the deuteron 5-wave 
only. The inclusion of the D-wave in the convolution with 1 S'o NN final state interaction is done 
like in the expression for diagram (al). Note that this diagram is antisymmetrized as well, but in 
this case the antisymmetrization reduces just to an additional factor of two, so we have written the 
resulting expression explicitly. The same is done for all diagrams with the half off-shell NN final state 
interaction, i.e. for the diagrams (b2), (c2), and (d2), and corresponding diagrams with pion higher 
partial waves. 

The properly anti-symmetrized amplitude Mnn for the S— and P-wave NN interaction with 
isospin 1 in the plane wave basis is 

(a' 1 a / 2 ,p'|m( j e;)|a 1 a 2 ,p) =2 ( x t, ^ x *,) m° nn (p,p',e) (xi; ^ XA 2 ) 

+2 (xi, <Ja^j= xl^j mU&p'^U (xl 2 ^jfp xx?) ,(23) 

where the first (second) term corresponds to the spin-singlet (spin-triplet) NN states. Here A' and 
A stand for the spin projections of final and initial nucleons, respectively, and xx for corresponding 
spinors. The amplitudes M.jy N in this formula are related to the standard NN partial wave amplitudes 
fg L , where J denotes the total angular momentum, L the angular momentum in the initial and final 
state (we do not consider couplings between partial waves with different angular momenta, e.g., the 
coupling of the 3 P 2 to the 3 F 2 partial wave), and S the total spin, through 

M° NN (p,p',E) = 16irM N f°°(p',p,E) 

MU&PiE)^ = 167rM^(2J + l)7>^//V,P,£) • (24) 

J 

The fg L (p',p,E) are related on the mass shell to the partial phase shifts via 

fi L (p,P,p 2 /M N ) = ±- (exp(2^ L ) - 1) . 
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The projection operators for the various total angular momenta with L = L' = 1 and for the S = 1 
read 

K/3 = PaPP , 

K P = \ e aXT ^ PT p'xP P = \ (5 af s(p ■ P') ~ P'pPa) , 
= ^ Uapip ■ P') + PpPa ~ ^PaPp^j , 

where p = p/\p | and analogously for p'. These projectors are normalized such that 



/ 



and the total neutron-neutron cross-section reads then 

<t = \-^Y. (2J + l)\f J s L (p',P,E)\ 2 • (25) 



2 JLS 



Note that the factor of one half occurs in this expression because the nucleons are identical. 



3. Diagram (bl). 

The real part of the diagram (bl) renormalizes the bare vertex in the leading diagram (al) resulting 
in the experimentally observed value for ^p — > 7r + n process. Thus, the real part of this diagram 
is already included in the expression for diagram (al). The imaginary part of (bl) is 

M bl = iC L0 87r ^ 1+ ™ W ^ MN ^ 2 u{p^/2) k 7rNlX \(ve~ ( )(ae D )^x*2 ~ (1 <-> 2), (26) 

where k^N^lM^k^— m 7r pj|/(Mjv+m vr ) is the magnitude of the relative momentum of the final pion 
and the final nucleon with momentum pi and Clo = ^V^NdnN^/ (v^)- 



4. Diagram (cl). 



where the integral I C1 is 



/ _fp_ ujp+q^/2) 
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5. Diagram (b2) 



M b2 = C LO gj^p^jeD) (x\ ^= X2) . (29) 

with 

i f SpdH M° NN (p+k n /2,p 12 , E l2 ) ujp+q^/2) 

The integral with the ttN loop, i.e. over d 3 l, is divergent and has to be renormalized (cf. discussion 
of diagram (bl)). After renormalization it takes the form 

where m^N = m^M^ / '(MN+m n ), = MN(MN+m 7T ) / '(2M]\r+m 7T ), and 

K(Q,p) = y/2m nN {Q- P y2 m ) . 

Note that for negative arguments the square root needs to be replaced by i \/2m 7T ^(p 2 /2fi]y—Q). Thus, 
we get a purely real, non-vanishing contribution from diagram (b2) even at production threshold. 



6. Diagram (c2) 



M C2 = C LO M^ptcA^D) (x\ ^= X*2 ) , (32) 



with 



The integral I C2 develops an angle dependent three-body singularity as soon as we move away from the 
production threshold. Fortunately, for non-relativistic pions, it transformed into an angle independent 
one. Indeed, rewriting the denominator in I C2 that creates the three-body singularity, as follows 

p 2 I 2 {P+1?_„ P 2 (l+^P) 2 



Q-^T7 ^77 = Q- 



2M N 2M N 2m n 2/i N 2m nN 

and shifting the integration variable / by p, we get the integral over d 3 l in Eq. (33) as 



/ 



dH u(f+&/2) f dH u(l + g 7 /2 -Jg*g 

- 2m 7rN / — — " , (34) 



(2vr) 3 __^ pL_M! + ,-n J ( 27 3 ^ 2 (Q,p)-/ 2 + iO ' 

2M*, 2Mk, 2m, * U 



2M N 2M N 2m T 

from which one can immediately see that the moving three-body singularity turned to the frozen one 
- the pole position does not depend on I angles. This is possible, because in our reaction the initial 
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two nucleons are always off the mass shell - they form a bound state. Therefore, there is no additional 
propagator which can cause a singularity. 

Furthermore, for non-relativistic pions and specific analytical parameterizations of the deuteron 
wave functions, the three body singuilarity can be integrated analytically. In particular, using the 
parameterization for the CD-Bonn deuteron wave function [32] (cf. appendix B) we may write 

f dH u{l + g 7 /2) 2m w 1 ( fa - iK(Q,p) + ir \ 

J (2vr) 3 q _ P _ p 2 {i+P) 2 , , n " 1 + mjM N ^ j 8mr ° g \fa - iK(Q,p) - ir J ' 

^ 2M N 2M N 2m, <U J 

(35) 

where r = i+m 1 /m n P ~ The integration in Eq. (35) has been performed similarly to the case 

of diagram (a2): we applied an inverse Fourier transform separately to the deuteron wave function 
and to the Green's function. Then the integration over dH gives a delta-function, which kills one 
of the remaining three-dimensional integrations and the integral reduces to a sum of integrals of 
the type considered, e.g., in Ref.[17]. The three-body singularities of the diagrams (dl) and (d2) 
can be handled analogously. Also for the NN final state interaction the integration over one of the 
loops can be performed analytically after applying the procedure outlined in Refs. [33, 34] to the 
NN amplitudes, for they can then be represented in the same analytical form as the deuteron wave 
functions, cf. appendix B for details. Moreover, it turns out that in diagram (d2) one can perform 
analytically two of the three loop integrations, thus reducing the nine-dimensional integral to the 
three-dimensional one. 



7. Diagram (dl) 



M * = c - ssb'*^ (*! 71 xi ) " (1 " 2) ' 



(36) 



with 



f d s pd 3 s M° NN (p+s/2, p 2 +s/2, E NN ) u{p+q/2) 
^ J (2„T M N E NN -(p 2 + s/2)^0 Q-^-^-lgW ( } 

where E^n = Q—s 2 /^^ with [i n = 2M^m i: / '(2M^+m 7r ). 



8. Diagram (d2) 



with 

h 2 



dH d 3 s d 3 p M° NN (p+k 7r /2,p 12 , E 12 ) 



M N J (2vr) 9 pf 2 -(P+^r/ 2 ) 2 +*° 
M° NN (l+s/2, p+s/2, E NN ) u{l+q/2) ^ 



M N E NN -(p+s/2) 2 +iO n j2 s2 (f+g) 2 , - n 

^ 2M N 2m, 2Mjv " r ' u 
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9. Corrections to diagrams (al) and (a2) from higher pion partial waves. 



We start from the explicit expression for the pion p-wave contribution (below labeled as (n—p)) 
stemming from diagram (b) of Fig. 4. To implement this, the term (<re 7 ) in the expressions for diagram 
(al) and (a2) of Fig. 1 (see above) needs to be replaced by ((?g)(e* 7 A: vr )/m vr , where q = <f 7 /|cf 7 |. Thus 
we get 

Mir P) =CLO —{hK) X \{aq){ae D )^ X l u{q 2 ) - (1 - 2) 
for the diagram without final state interaction and 

aA*-p) n 1 c t \ 1 [ d3p M nn(p+^/'Z,Pi 2 ,E 12 ) / t o 2 * 

for the diagram with NN final state interaction. 

Note, that when the pion is in a p-wave, only those terms where the NN final state is in an S— 
wave are to be considered. The simultaneous appearance of two p-waves in the final state is strongly 
suppressed by the centrifugal barrier. Our numerical calculations confirm this statement. 

Using the vertex V^j NN given by Eq. (16) as input for the diagrams (al) and (a2) of Fig. 1 one can 
get the corresponding contribution from the s- and u-channel nucleon pole diagrams (cf. diagrams (c) 
in Fig. 4) to our reaction as follows 

Mnuc ~ su = ~2§~k; { {Anuc ~ su + pnuc ~ SU) ( xti 72 xl ) + §nuc ' su ( x * d 72 xl ) } ' (40) 

where 

A nuc-su = u ^ ^2{e^p 2 )(Ke D )+^ p + fi n ) ((? 7 4)(9 7 eD)-(? 7 4)(e 7 ?D)) } +(1 ~ 2), (41) 



gnuc su = u (^ 2 )|2(e 7 p 2 )(4xe D ) + (/x p -/x n )(4(e 7 xg 7 ))eD 

+ (fJ-p + fJ-n) ((e 7 4-)(g 7 x e D ) - (q 7 k n )(€y x e D fj } - (1 <-> 2), (42) 
and F nuc ~ su stands for the contributions from the diagrams with NN final state interaction: 



pnuc su = ^ + ^ ^^fafryfakn)^^ F 1 - (k n € 1 )(k n e D )F 2 , 

where 



Fi = — ( u&rt,l% 

2M N J (2vr)3 p 2 2 _^ +k7T / 2 y+iO 



p 2 = 1 I d 3 p M% N (p+kj2,p 12 ,E 12 ) ^ i ^, 2 ^ i {p + K/2){K 12-^/2) 



2M N J (2vr)3 p j 2 - (^+4/2)2+^0 " \ (4/2 - «f 7 /2) 2 
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B The NN wave functions 



For the deuteron wave function, and for the nn scattering amplitudes that appear in the final-state 
interaction, we take those of the (charge dependend) CD-Bonn NN potential [32]. In particular, we 
utilize the analytic parameterization of the deuteron wave function provided in Ref. [32] which is 
given by 

u( P ) = ^Y. c i/(p 2 + m f ); w (p) = ^Y. D oi^p 2 + ( 43 ) 

with parameters listed in Table 20 of that reference. The wave function is normalized according to 

(1 " P (u 2 (p)+w 2 (p)) = l. (44) 



(2vr) 2 



With this parameterization some of the diagrams can be evaluated analytically. In order to facilitate 
also an analytic evaluation of the diagrams involving the nn scattering amplitude the CD-Bonn poten- 
tial in the relevant partial waves ( 1 S'o, 3 Po, 3 P\, 3 P2) is cast into a separable representation by means of 
the so-called EST method [34]. The resulting rank-1 separable interactions exactly reproduce the on- 
and off-shell properties of the CD-Bonn potential at the chosen approximation energies {El^ = MeV 
for 1 So and E^ = 30 MeV for the P waves) [34] and they provide also an excellent approximation in 
a broad neighborhood of these energies. The form factors g(p) of these separable representations, that 
consist of the scattering solutions of the CD-Bonn potential at the specified approximation energies 
[34], are parameterized in analytical form, 

s(p)=;>>/^ 2 +& 2 )' ( 45 ) 

i 

for 1 Sq and 

g(p)=Y,dP/(P 2 + fi) 2 ' (46) 

i 

for the P waves and the scattering amplitude is then given by 

f(p,p';k)= ^ 2 Mn 9 (p)9(p>) 

Here the positive sign pertains to the 1 S'o, 3 Pq, and 3 P2 partial waves, and the negative sign to the 
3 Pi partial wave. The parameters q and (3i for each partial wave are listed in Table 3. 
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